Energy 72 (2014) 188-194 



ELSEVIER 


Contents lists available at ScienceDirect 

Energy 

journal homepage: www.elsevier.com/locate/energy 


BIE1BI 



Volatile species release during torrefaction of biomass and its 
macromolecular constituents: Part 2 - Modeling study 

Timothee Nocquet l , Capucine Dupont ' , Jean-Michel Commandre , 
Maguelone Grateau a , Sebastien Thiery , Sylvain Salvador 

a CEA Grenoble — LITEN — Laboratory of Technology of Biomass, 17 rue des Martyrs, 38054 Grenoble Cedex 9, France 
b CIRAD, Biomass Energy Unit, 73 avenue fean-Frangois Breton, 34398 Montpellier Cedex 5, France 
c MINES ALBI, RAPSODEE, FRE CNRS 3213, France 



CrossMark 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 20 September 2013 
Received in revised form 
4 April 2014 
Accepted 3 May 2014 
Available online 2 June 2014 


Keywords: 

Biomass 

Cellulose 

Hemicelluloses 

Lignin 

Model 

Torrefaction 


Based on the torrefaction experiments carried out on beechwood and on its macromolecular constituents 
- cellulose, lignin and xylan - in part 1 of this paper, a simple torrefaction model was developed to 
predict both solid and the eight main volatile species yields versus temperature, residence time and 
biomass macromolecular composition. This model consisted in the superposition of kinetic “sub-models” 
describing the torrefaction of each constituent. The interactions observed in part 1 between constituents 
were then taken into account by means of an empirical factor related to cellulose decomposition rate. 
The model reproduced quite well the experiments and enabled to show the influence of each constituent 
on torrefaction yields. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

The general context of the study is detailed in part 1 of this two- 
paper series. The reader is invited to consult the experimental part 
(part 1) before considering the following modeling part (part 2) to 
better understand the study. 

As mentioned in part 1, at the moment, torrefaction mecha¬ 
nisms are still poorly understood [1]. There are only few kinetic 
models dedicated to biomass torrefaction [2-4]. These models are 
mainly inspired from pyrolysis models, for which an abundant 
literature is available, as reviewed by Di Blasi [5]. Pyrolysis and 
torrefaction models both describe biomass thermal degradation, 
but in different ranges of temperature, i.e. the range of temperature 
in pyrolysis is higher than that associated to torrefaction - typically 
higher than 300 °C. Hence, the reactions and mechanisms involved 
are expected to be — at least partly — different and the substantial 
work done to develop pyrolysis models cannot be directly applied 
to torrefaction conditions. 
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As solid is the main product of torrefaction, the objective of 
torrefaction models is generally the description of solid mass loss 
versus time and temperature. Only Bates and Ghoniem have 
recently developed a model able to predict the composition of 
volatiles [3]. As highlighted by the authors, their model suffers from 
the lack of experimental data for validation, with only one study 
available about volatile species emission during torrefaction [6]. 
Moreover, although it has been shown that both solid yield and 
volatile species yields were highly influenced by the biomass type 
[6-8], the torrefaction models mentioned above generally do not 
describe the influence of biomass type. As most pyrolysis mod¬ 
els,they consider biomass as a global solid and are thus feedstock- 
specific. Only Rousset et al. [2] have developed a model which 
could be applied to various biomasses. Their approach was based 
on an approach previously developed in pyrolysis by Koufopanos 
et al. [9]. This approach consists in considering the summative 
contribution of the biomass macromolecular constituents, which 
are known to react differently under the effect of heat [10]. Un¬ 
fortunately, this model was limited to the description of solid mass 
loss versus time. Moreover, the model physical meaning was 
disputable, since its parameters were not derived from experi¬ 
ments on constituents. 















T. Nocquet et al. / Energy 72 (2014) 188-194 


189 


Nomenclature 

V! 

volatile species produced by step 1, - 



V 2 

volatile species produced by step 2, - 


initial solid, — 

V* 

complementary species, - 


intermediate solid produced by reaction 1, - 

y A 

Yield of remaining initial solid, wmf% 

^2 

torrefied solid produced by reaction 2, - 

VaI or A2 

yield of solid produced by step 1 or step 2, wmf% 

C S 

variation of mass percent of constituent s, wmf% 

Vs 

total yield of solid, wmf% 

Ea 

activation energy, J.moH 1 

vv 

total yield of volatiles, wmf% 

F 

function describing the difference between model and 

VVl or V2 

yield of volatiles produced by step 1 or step 2, wmf% 


experimental yields, wmf% 

V^sum i 

total yield of species i, wmf% 

k 

kinetic parameter, s -1 

«ivi or v 2 parameter describing the production of species i in 

ko 

preexponential factor, s 1 


step 1 or 2, — 

f^lor 2 

sum of kinetic parameters associated to step 1 or 2, s -1 

(Oj 

initial mass fraction of constituent j, wmf% 

R 

universal gas constant, J.mol _1 .K _1 

w% 

mass percent 

T 

temperature, K 

wmf% 

moisture-free mass percent 

t N 

interaction factor, - 




In part 1 of this paper, an original set of experimental data was 
supplied about solid mass loss and volatiles production during 
torrefaction of beechwood, of its macromolecular constituents, i.e. 
cellulose, xylan and lignin, and of their mixtures. These data 
showed the ability of the additive law applied on constituents to 
describe torrefaction up to 250 °C and the existence of interactions 
between constituents at temperatures of 280 °C and 300 °C, maybe 
due to the ramification of cellulose fragments by lignin and xylan 
radicals produced at the beginning of torrefaction. 

Based on the existing work in literature and on part 1 of this 
paper, the objective of part 2 is to derive from the experimental 
results a simple model able to predict both solid yield and the main 
volatile species yields during torrefaction versus temperature, 
residence time and biomass macromolecular composition through 
the use of a corrected additive law on biomass constituents. 

First, the model principle will be described. The parameters 
determination will then be detailed. To end, the model perfor¬ 
mances will be discussed as well as the influence of biomass 
macromolecular composition on the results. 


2. Model development 

2.1. Preambue 

The model assumptions rely on the large set of experimental 
results obtained in part 1. The model has therefore some physical 
basis and it may also be a useful tool to better understand the 
torrefaction processes. However, it does not claim to be an exten¬ 
sive phenomenological description of biomass torrefaction. Indeed, 
this semi-empirical model mainly aims at being a simple tool that 
may be integrated in reactor models for further process design. 

2.2. Model description 

The general principle of the model is described in the scheme of 
Fig. 1. The input data are the contents in biomass constituents, i.e. 
cellulose, lignin and hemicelluloses. Hemicelluloses are repre¬ 
sented by two constituents: xylan and acetyl groups. This choice 
enables to describe the production of acetic acid by the hydrolysis 



Fig. 1 . Principle of the model. 
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Table 1 

Kinetic parameters value for the different constituents. 


Constituent 

Kinetic parameter 

ko (s- 1 ) 

Ea (kj.mol 1 ) 

Cellulose 

l<A\ 

1.32 x 10 12 

171 


kv i 

1.06 x 10 21 

268 


/<A2 

2.10 x 10 33 

392 


/<V2 

3.09 x 10 18 

230 

Lignin 

/<A1 

2.33 x 10 3 

92 


ky i 

1.36 x 10 5 

88 


/<A2 

1.03 

38 


kv.2 

3.41 x 10 4 

93 

Xylan 

/<A1 

1.18 x 10 7 

96 


fcvi 

1.67 x 10 10 

133 


/<A2 

7.21 x 10 8 

127 


kv.2 

2.50 x 10 9 

136 

Acetyl groups 

/<A1 

1.61 x 10 2 

59 


fcvi 

6.46 x 10 6 

110 


/<A2 

5.32 x 10" 2 

9 


kv 2 

4.22 x 10 12 

149 



Fig. 2. Comparison between experimental solid mass loss and model prediction at the 
four temperatures tested. 


of acetyl groups, which are absent from the extracted xylan, as 
explained in part 1. Besides, as explained also in part 1, ash and 
extractives are not taken into account. 

For each constituent, a kinetic sub-model was developed. As 
most existing torrefaction models, each sub-model was based on 
the semi-global mechanism proposed by Di Blasi and Lanzetta [11 ]. 
It enables to describe the two-phase process observed in the evo¬ 
lution of mass loss versus time during torrefaction. Indeed, this 
mechanism consists in a series of two sequential reactions, as 
shown in Equation (1 ). First, two parallel reactions occur. One gives 
rise to volatile species called V\, the other one gives rise to an in¬ 
termediate solid A\. This intermediate solid is then assumed to 
react, as previously, through two parallel reactions. One gives rise 
to volatile species named V 2 and the other one to the torrefied solid 
A 2 . Each reaction x was assumed to obey a first order kinetics and 
the associated kinetic parameters k x were assumed to follow the 
Arrhenius law (see Equation (2)). 



where k 0x is the preexponential factor, Ea x the activation energy, R 
the universal gas constant and T the temperature. 

Based on this mechanism, the solid mass Y s is equal at each 
moment to the sum of the A, A\ and A 2 masses Ya, Ym and Ya 2 
(Equation (3)) and the total volatile yield is Yy equal to the sum of V\ 
and V 2 yields Y V i and Y V 2 (Equation (4)). 


y s = y A + Ym + ^2 

(3) 

Vy = Vyi + y V2 

(4) 


As reactions are assumed to obey first-order kinetics, the yields 
in solid Y s and in volatiles Y V i and Y V 2 are given by Equations 
(5)—(7). The obtention of these equations is not mentioned here as 
it can be found in the papers of Di Blasi and Lanzetta [11] and of 
Prins et al. [12]. 




(K Al^l - 

V ) 


exp(-/C 1 t) 


-/<A1 K 2 + /<A1 /<A2 


K 2 (K 2 -K i) 


exp(-/C 2 t) + 


f<Al f<A2 

K,K 2 


(5) 


y V i =^ L(i - exp( -Xit )) (6) 

Yv2 = KikJk ?- l< 2 ) (K] - I<2 - K ' ex P(- K 2 [ ) + 1 <2 exp(-/Ci t)) 

(7) 

with: 

K i — k ai + /<vi (8) 


I<2 = k A2 + /<V2 ( 9 ) 

Volatile species from the two reactions were assumed to be 
composed of the eight main species produced during torrefaction 
experiments: carbon dioxide, carbon monoxide, water, formalde¬ 
hyde, acetic acid, furfural, methanol and formic acid. To describe 
the volatiles composition, a parameter cq x was associated to each 
species i for each reaction. The yield of the species was then given 
by Equation (3). 

Vsum i — a iVi Vvt + a iV 2 Vy 2 (10) 

where Y sum ;, Y V i and Y v2 are the global yield of species i, the yield of 
volatiles by reaction 1 and the yield of volatiles by reaction 2 
respectively and with 


Table 2 

Value of the coefficients associated to the production of volatile species for each constituent. 


Constituent 

Carbon dioxide 

Carbon monoxide 

Water 

Formaldehyde 

Furfural 

Formic acid 

Methanol 

V x 

Cellulose 

0.028 

0.008 

0.426 

0.120 

0.001 

— 

— 

0.417 

Lignin 

0.045 

0.020 

0.350 

0.470 

0.004 

- 

0.027 

0.084 

Xylan 

0.233 

0.077 

0.430 

0.025 

- 

0.023 

0.042 

0.170 
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Table 3 

Comparison between experimental and model yields at 220, 250, 280 and 300 °C for the main volatile species. 


Temperature (°C) 

Yield (wmf%) 

Carbon dioxide 

Carbon monoxide 

Water 

Formaldehyde 

Furfural 

Formic acid 

Acetic acid 

Methanol 

220 

Experimental 

1.48 

0.50 

3.42 

1.35 

0.01 

0.13 

0.90 

0.31 


Model 

0.67 

0.10 

2.22 

1.59 

0.00 

0.00 

0.95 

0.47 


Absolute error 

0.81 

0.40 

1.20 

0.24 

0.01 

0.13 

0.05 

0.15 

250 

Experimental 

2.93 

0.99 

7.29 

2.58 

0.02 

0.26 

3.37 

0.60 


Model 

2.43 

0.38 

6.21 

3.46 

0.05 

0.77 

3.18 

0.77 


Absolute error 

0.51 

0.61 

1.09 

0.88 

0.04 

0.51 

0.19 

0.17 

280 

Experimental 

4.20 

1.40 

16.97 

5.97 

0.05 

0.33 

5.45 

0.79 


Model 

3.43 

1.11 

12.45 

5.84 

0.17 

1.53 

5.32 

1.09 


Absolute error 

0.77 

0.29 

4.52 

0.13 

0.13 

1.20 

0.13 

0.30 

300 

Experimental 

4.93 

1.64 

23.24 

8.31 

0.07 

0.36 

5.90 

0.90 


Model 

4.72 

1.82 

21.01 

8.15 

0.33 

2.04 

6.14 

1.50 


Absolute error 

0.21 

0.19 

2.23 

0.16 

0.27 

1.68 

0.25 

0.60 


E = E aN 2 = i o i) 

i i 

As the composition of volatiles was found to be roughly the 
same under the range of temperatures investigated during exper¬ 
iments, these parameters were assumed to be constant whatever 
the temperature. Besides, as a first approximation, the volatiles 
composition was assumed to be constant over time and the pa¬ 
rameters am and am were thus assumed to be equal. Since it was 
observed that each constituent did not produce all species, some 
constants were taken equal to zero: 

• Cellulose only gives rise to carbon dioxide, carbon monoxide, 
water, formaldehyde, furfural; 

• Lignin only produces carbon dioxide, carbon monoxide, water, 
formaldehyde, furfural, methanol; 

• Xylan only produces carbon dioxide, carbon monoxide, water, 
formaldehyde, methanol, formic acid; 

• As explained above, acetyl groups produce acetic acid. 

As all volatile species were not measured in experiments, a 
complementary species V x was added to close the mass balance. 

To obtain the biomass torrefaction model, the additive law was 
applied by summing the contribution of each constituent sub¬ 
model. As discussed in part 1 of this paper, interactions could be 
observed between cellulose and the two other constituents at 
280 °C and 300 °C, maybe due to the ramification of cellulose 
fragments by lignin and xylan radicals produced at the beginning of 
torrefaction. To describe these interactions, a corrective interaction 
factor tN was added to the kinetic parameters k ai, /<A 2 , /<vi and km in 
cellulose sub-model. 

2.3. Parameters determination 

To determine the values of model parameters, a function F was 
defined. This function F characterizes the difference between 
experimental yield and predicted yield (see Equation (12)). F was 
minimized for the four temperatures investigated during the 
isothermal period of the experiments. 

F = £| modelexp| (12) 

t 

23A. Sub-models parameters 

To determine for each sub-model the kinetic parameters, i.e. the 
preexponential factor 1<o x and the activation energy Ea x , the func¬ 
tion F was defined as the absolute value of the difference between 
solid yield versus time obtained in the thermogravimetric experi¬ 
ments (see part 1 of this paper) and solid yield calculated from 
Equation (5). 



Fig. 3. Model predictions of solid mass loss versus time at 220 °C, 250 °C, 280 °C and 
300 °C for different contents in a) cellulose; b) lignin and c) hemicelluloses. 
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Fig. 4. Model predictions of volatile species yields for different contents in cellulose after 3 h of torrefaction at a) 220 °C; b) 250 °C; c) 280 °C; d) 300 °C. 


To determine for each sub-model the parameters describing 
volatiles composition, F was defined as the absolute value of the 
difference between total yield of volatile species i obtained in ex¬ 
periments inTORNADE (see part 1 of this paper) and that calculated 
from Equations (6), (7) and (10). 

The values obtained are given in Table 1 for the kinetic param¬ 
eters and in Table 2 for the constants describing volatiles 
composition. 

It appears hazardous to make a direct comparison between the 
kinetic parameters found in this study and the kinetic parameters 
of pyrolysis models. Indeed, as explained in introduction, these 
models are valid in a range of temperature different from those of 
torrefaction and may therefore describe other phenomena. Simi¬ 
larly, it seems difficult to compare kinetic parameters associated to 
different mechanisms, i.e. one-step or two-step. However, a com¬ 
parison could be made with the kinetic parameters found by 
Repellin et al. on hemicelluloses with the same two-step mecha¬ 
nism at the same temperature level [4]. The activation energies 
found by Repellin et al. were 30-80 kj moH 1 lower than the pre¬ 
sent ones and surprisingly the preexponential factor he found for 
the second reaction of volatiles production was equal to zero. Such 
differences must be due to the compensation effect occurring in the 
kinetic parameters determination: different couples of pre¬ 
exponential factor and energy activation may fulfill the numerical 
criteria. 

2.3.2. Factor describing the interactions between constituents 

To determine the interaction factor tn, F was defined as the 
absolute difference between solid yield versus time measured 
during mixtures tests and solid yield predicted based on the 
behavior of constituents alone. 


The value was found to be nearly the same when obtained from 
tests on cellulose—xylan or on cellulose—lignin. This shows that the 
interactions have the same effect on the mixture mass loss 
evolution. 

This value varied with temperature: it was of 1 at 250 °C, of 
0.54 at 280 °C and of 0.32 at 300 °C. This means that the in¬ 
teractions effect tends to become stronger with temperature with 
a more significant decrease in cellulose decomposition rate, that 
is, approximately of a factor 2 at 280 °C and of a factor 3 at 
300 °C. 

3. Results 

3.1. Model performances 

As shown in Fig. 2, the model was able to predict with satis¬ 
factory accuracy the solid mass loss versus time within the range of 
temperatures tested. Indeed, the absolute difference between ex¬ 
periments and model was always below 3%. 

The experimental and model yields are compared for the main 
volatile species in Table 3. The differences between experiments 
and model were relatively similar for all temperatures but varied 
according to species: the relative difference was of 35% or less for 
carbon dioxide, water, formaldehyde, acetic acid and methanol. It 
was slightly higher (50%) for carbon monoxide. Hence, for all these 
species, the prediction could be considered as quite good as it en¬ 
ables to get a correct estimation of the amount of products. The 
results were not so satisfactory for two minor species: formic acid, 
with a difference of 200% and furfural, with a difference of 300%. 
However, such high difference may be explained by the very low 
amounts of these species in the experiments. 
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3.2. Influence of the initial biomass composition 

Experiments have not been performed on other biomass sam¬ 
ples with different macromolecular composition yet. The model 
was then used to investigate theoretically the influence of biomass 
macromolecular composition on the torrefaction process. Model 
calculations were performed by varying the amounts of cellulose, 
hemicelluloses and lignin in the initial biomass. The values were 
chosen in the limit of the range typically encountered in literature 
about biomass 13]: 

• Cellulose from 35 to 60 w% 

• Lignin from 15 to 35 w% 

• Hemicelluloses from 15 to 35 w%. 

Note that these three parameters are not independent. It was 
chosen to compensate the variation of a constituent by the 
variation of the two other ones, the ratio between them 
remaining unchanged. Hence a decrease of the constituent 1 of x 
w% resulted in an increase of the constituents 2 and 3 of C 2 X w% 
and c 3 x w% respectively, with C 2 and c 3 given by Equations (13) 
and (14). 


to a higher production of water when cellulose content 

increases. 

4. Conclusion 

Based on the experimental results obtained in part 1 of this two- 
paper series, a semi-empirical model was developed. This model is 
the first literature attempt to predict both solid and the main vol¬ 
atile species yields (carbon dioxide, carbon monoxide, water, 
formaldehyde, acetic acid, furfural, methanol and formic acid) 
produced during biomass torrefaction versus temperature, resi¬ 
dence time and biomass macromolecular constituents. The addi¬ 
tion of kinetic sub-models developed for each constituent - 
cellulose, lignin, xylan and acetyl groups - gave satisfactory results 
provided a corrective factor was added on cellulose decomposition 
rate at temperatures above 280 °C. This corrective factor took into 
account the interactions between cellulose and the two other main 
constituents. In the future, this model may be integrated in reactor 
model. 
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^2 


Cl>2 

Gl> 2 + w 3 
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c 3 


Gl > 3 

( jl >2 + W 3 
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where coj is the initial mass fraction of the constituent j in biomass. 
As illustrated in Figs. 3-6, the main results were: 

• An increase of lignin content has no effect on solid yield; it 
mainly leads to an increase of formaldehyde yield. 

• An increase of hemicelluloses content impacts the mass loss that 
becomes higher. This increase is logically visible at the begin¬ 
ning of torrefaction, when most hemicelluloses are degraded. 
This gives rise to an increase of carbon dioxide, carbon mon¬ 
oxide and water. 

• A higher cellulose amount leads to a decrease of solid mass loss 
and volatiles yield at 220 °C and 250 °C, which is logical since 
this constituent is not attacked at these temperatures. 

• The mass loss of solid with the highest cellulose amount is 
higher at 280 °C and 300 °C after 60 min of torrefaction than 
the mass loss of solid with the smallest amount. This is related 
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